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digital simulation. This model is combined with that of a mariner hull. 
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I. INTRODUCTION 


The propulsion system of a ship is a complex, nonlinear system with 
several internal feedback loops. The performance of the ship itself is 
dependent, of course, on the capabilities and limitations of the pro- 
pulsion plant, which may limit overall ship performance so that the 
.capabilities of the hull.design are not realized. 

The first objective of this thesis is to model, in detail, a 
typical ship propulsion system, and to simulate it in the digital 
EL After validation of the computer model a simulation of the 
Mariner hull is to be coupled to the propulsion plant so that control 
system studies may be undertaken. As part of the validation procedure, 
behavior of the shaft speed was investigated with and without the 
governor. 

The second objective was Ee. the computer model as part of a ship 
feedback control system. For this purpose it was decided to explore the 


use of a feedback loop to maintain constant ship speed in turns. 


Liam 360/67 DSL 360 LANGUAGE at W. R. Church Computer Center, Naval 


Postgraduate School. 


p 





II. PROPULSION PLANT DYNAMICS 


A. SHIP PROPULSION EQUATIONS 

The propeller angular speed is denoted by n. 

The basic equation describing the angular acceleration of a propul- 
sion shafting system is: 

"EQ = 2п14п/4г 


Where: XQ = Summation of torque = Qe-Qp 


I = Polar mass moment of inertia 

n = Propeller angular speed 

Qe = Shafting torque 

Qp = Propeller reaction torque 
then dn/dt = ди EQ 

211 


For the typical gas turbine, the experimental data on the power 
turbine torque Qe versus the power turbine speed is shown in Fig. Al. 
The torque map of Fig. Al allows the engine torque to be determined if 
the engine speed and fuel flow rate Wf are know. However, this torque 
representation is correct only for steady state contitions. 

Propeller speed n(rps) is related to the engine speed N3(rpm) 


through the reduction gear ratio Kg 


N3 = 60Kg.n 


The block diagram for the ship propulsion system is given in Fig. A2. 
The power turbine speed can be controlled by a governor acting on 


the gas flow. The changes in fuel flow and speed that occur can be 
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assumed small enough to permit the use of constant coefficients in the 
dynamic equations. 
Using the time lag of the governor,t, the block diagram for the 


Governor is shown in Fig. A3, 


Where: n* Rotational speed order 


n = Rotational speed 
ner = n* —n = Speed error signal 
Kf = Fuel flow rate/RPM 
T = Time delay 
W£ = Fuel flow rate 


GOVERNOR + КЕ/(1+15) 


Then the propulsion plant is shown in Fig. A4. 
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Figure A2. Block Diagram of Ship Propulsion System 
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Figure 


A3. 
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Block Diagram for the Governor. 
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Figure A4. The Propulsion Plant with Governor 
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III. HULL EQUATION AND SIMULATION 


WAKE FRACTION: W 
The difference between the ship speed and the speed of advance of 
the propeller Vp is called Wake speed. 
The Wake fraction is defined as 
W = (V-Vp)/V 
thus Vp-V(1-W) 
THRUST COEFFICIENT: Ct 
TORQUE COEFFICIENT: Cq 


For the propeller in open water 


ЕЕЕ === (1) 


O A (2) 
пон) 


Wnere r = Thrust 
Q = Torque 
p = Water mass density 


Ct = Thrust coefficient 


Cq Torque coefficient 


а 
n 


Propeller diameter 


From (1) and (2) 


T 29 (7007) 


Sa gale Cie в) 


© 
! 
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The experimental data on the Wake fraction versus the speed of the 
ship is shown in Fig. Bl. 

The Thrust coefficient Ct and Torque coefficient Cq versus the 
Second modified advance coefficient are shown in Fig. B2 and Fig. B3. 


Where Vp = V(1-W) 


And o DVDS pA) 


Total-ship Resistance, Rt, is a computer look-up table giving Rt 


versus Ship speed V (fig. B4). 
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Figure B2. Thrusting Coefficient Ct versus Second 


Modified Advance Coefficient o 
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Figure 84. Ship Resistance versus Ship Speed V. 





IV. COMBINED PROPULSION PLANT AND HULL EQUATION 


A. POWER TURBINE WITHOUT GOVERNOR: 


Nonlinear functions such as Wake fraction, Thrust deduction fraction, 


Thrust coefficient, torque coefficient, ship resistance and engine torque 


were stored in look-up tables during the computer solution. 


Using the DSL/360 language these equations were solved by means of 


an integration and updating from table look-up every 0.1 sec until the 


steady state was reached. 


Ere. 
Fig. 
Fig. 
Fig. 
Fig. 


Fig. 


B5. shows 
B6. shows 
B7. shows 
B8. shows 
B9. shows 


B10.shows 


propeller shaft speed (rpm) versus time 
ship speed V versus time 

propeller action torque versus time 
shaft torque Qe versus time 

Cq versus time 


Ct versus time 


For a step change in fuel flow rate 


Wf = 7000. + 3000. step (0.0) 


The shaft torque decreases and goes to steady state. The propeller 


action torque increases, overshoots at 10 sec. then goes to steady state. 


The propeller shaft speed increases, goes to steady state in about 80 sec. 
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Figure Bl10.- Torque Coefficient Cq versus Time 
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B. POWER TURBINE GOVERNING: 

Turbine speeds with governor are shown in Fig. A4 and compared with 
the responses that occur without a governor. 

For small step changes in power demand at 7000 Lb/hr fuel flow rate, 
the governor, which has a rate of 40 Lb/hr. RPM and time lag 1/; =1/5=0. 2 
sec., causes the power turbine response to be reduced from 80 sec. to 


20 sec. with an overshoot of 3.92% of the final change in speed. 


RPM command 

(52120512 ӘТЕРІО;) 
Now changing 
= 45. RPM (Fig. Cl) 
= 40. RPM (Fig. C2) 
30. RPM (Fig. C3) 


= 20. RPM (Fig. C4) 


a E 
|| 


= 10. RPM (Fig. C5) 


Percent Overshoot as a function of rudder angle is shown on Table l. 
The block diagram for the combined propulsion plant and hull 


equations were shown in Appendix B. 
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V. EQUATIONS OF MOTION OF THE HULL 


A. NEWTON'S LAW OF MOTION: 

Newton's Law of motion for a rigid body can be written as two 
equations: 

FORCE EQUATION 


— —,— ы 
F = Force on body = d(Momentum)/dt 


4 (ш) 


EE 


Where Че Velocity of body 


MOMENT EQUATION 


--> 
M = Moment acting on a body 


d(Angular momentum) 
dt 


а(19) 


ас 





F Where I = Moment of inertia 


52 = Angular velocity 


В. SIX DEGREES OF FREEDOM: 


А, Ne | V % 


Ax 


X Jy 
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- 


Redrawing: 





ма 


€ OF y 


2 


Vertical axis, positive downward 
Tranverse axis, positive staboard 
Longitudinal axis, positive forward 


Roll angle 

Positive rotations are indicated 
Pitch angle 

in the sketch above 
Yaw angle 








== — 
F = Er 3 mU + X(mU) 
=> ur > > — 
{ = а(н) = Н+. Ох Н 
dt 


m: Mass of ship 


— 
U 7 uitvjtwk (3) 
u: Rate of surging 


v: Rate of swaying 


w: Rate of heaving 
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——À5 
R=pi+tqj+rk (4) 


R: Angular velocity 
p = Фф = Rate of roll 


O - Rate of yaw 


а 
|| 


r = Ф = Rate of pitch 


T = Xi + j + Zk (5) 


X: Hydrodynamic force along X axis 
Y: Hydrodynamic force along Y axis 


Z: Hydrodynamic force along Z axis 


=| 
|| 


Moment vector acting on body 


Ki Tp MJ + Nk (6) 
K: Rolling moment about X axis 
M: Pitching moment about Y axis 


N: Yawing moment about Z axis, 


Substituting equations (3),(4),(5), and (6) into equation (2) 
gives 


X = m(u + qw - гу - x. (a° + 20) 4 Y. (pq _ r) + Zo (Pr + q)) 
° 2 2 5 . 
^ Y =m(v + ru - pw - Y. (p + r”) + Z. (qr - p) + X. (qp + r)) 


. 2 2 . 2287 
2 m(w + pv - qu - Z. (P +q) + X.(rp - q) + 7 (та + р)) 


K = Тур + (1; - Ly) qr - m(Y (о + pv -qu) - Z. (v + ru - pw)) 
M = Lyq + (1, - І,) үр - m(Z¿(u + qw — rv) - X. (w + ру - qu)) 


М = is + (1, - Ту) pq - n(X.(v + r = pw) - Ү (и + ам - ту)) 


For the surface ship in calm water, roll, pitching and heaving are 


taken as zero. That means p=q=w=0 and most ships have Yo70. 
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Then equation (7) becomes: 


X - m(ü-rv-X T 

G 
Y = m(vtrutX,r ) (8) 
N = I,r-m(X,(vrur)) 


С. LINEARIZATION OF THE EQUATIONS OF MOTION: 


l. X Equation 


The nonlinear TAYLOR expansion for the X equation is 


Xz X* e [Xu +XyV+X2 + 5] 
2 
fs A (A3) 4 Хуу V Хоу vt Xy 


мы... ‚+ 2Х„®6 | 


(5) 


% 5 5 
+ 6 | аы СА Е Xow + a 
2 


5 Хау ОЎУ 3 Хан (40 + 6 Хр 0 


t6Xu 09 4 sls. " 


Note: 


1. The TAYLOR expansion of the equation above includes terms up 
to third order, since terms higher than third order are not significant. 
2. Since X is an EVEN function of v, of , of 
X(v) = 2 + ануч+ ав У... 
G, te by ДЕ + + ° 


24 сад Ste © 


Х(г) 


х(5) 


Then the odd powers of v,r,d are zero i.e. all Х У.Х х,Х%4 аке equal 


Zero. 
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Thus equation (5) becomes: 


| 2 2 2 
Х-Х- E Xuy (4) + E Хуу ЈЕ Хуа? +259 


+ Хр, (609 + X a (Au) 2 4 Xy 5 о + Ху X vj + X, 75 (6) 


Some of these terms are equal to zero (see Table 2) 


then X-X970 


2. Y Equation 


The nonlinear TAYLOR expansion for the Y equation 15: 


ЕУ + аи Фу + ул +48] 
Е Е ам A е + Yun e 
+2 У, |5 
4 € Б (44) + Y t Ха + Ууд” 
23 Yyy (DV +? Y + в t 


US ту A a | (7) 


NOTE: 


1. The terms higher than third order are not significant. 


2. Y is an odd function of v,r,ó 
EOS A р аб. 
тэтгэн 
е 5 
YO) - 446 pr 
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Then equation (7) becomes: 


~ 


DES ДИ Е Ж л. +)5д + 
2 2 3 5 
ve ee (Au) + > ју У + = Yann + L 465 
2 1. 2 22 2. 
Ша Да (Au) V ЩЕ Yıur (4407 t 9 4 +3 Yaw™ ty Myo? 


2001 2 
+š Wss YÍ 52 345 лд + Ум YZS 


Where the indicated terms are equal to zero (see Table 3) 


У-у + VA A Low P уа y na 
(8) 
3. N Equation 
The nonlinear TAYLOR expansion for the N equation is: 
М=№+ Мдм + Ми + Nat +N5d+ 
+5 | pu (OY Nw Ма + рад + 
12 N +... фо № | 
Nunu ли“ ЕМ + 177772: UP TM 
42 Nu (9 +35 Ма (АША... > Малта? 


+ 6 Muy, 38.7 + У 12201724 


Ци 
u 6 


(9) 


NOTE: 


1. The terms higher than third order are not significant. 


2. N is an odd function of v,r, £ 


8:7 + 95 > да! + . 


N(v) = 
N(r) = ОО МТ үр. + ER 
м8) = Кб КР КАС" 
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Then equation (9) becomes: 
N= Ny Au + МУ+ М + Мд + 
I 2 5, 2 
6 1797 (ди) + = Муу V + >> Мара, 1 - с Муза ЛБ 
2. 2 
+g Nuw (Au) V+ М№ах 207 ho... Мад 
The indicated terms are equal to zero (see Table 4) 


Then 


2 2 2 
N = Nv V rN; +M ó + z Мууу У?+ = Мад 7 Мур, УФ Ta Мами 2У 
(10) 


Combining the third order TAYLOR expansion for X,Y, and N with the 


dynamic response terms of the X,y, and N equation. 


|| 
O 


X equation X 


Y equation 


(m-Yv)vt(mX,-Y;) r = Y(u,v,r,8) 
N equation (1) 
(mX,-N-)vt(I,-N:)r 2 Ки м,т.8) 
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Using CRAMER'S rule to solve these equations 


ще 
|| 


Eu vir) 
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TABLE 2. ASSESSMENT OF THE COEFFICIENTS IN 
THE X EQUATION 


Variable Coefficient Series 60/Model 5.1.1 


m-Xu 0.1750 
Xu 

1/2 Xuu 
1/6 

1/2 

1/2 

1/2 

1/2 

Ша 

1/2 Х%ч 
Xvr+m 
Худ 

Xr 

Xvru 

Ху u 

Xr u 

Xo 


A 
U ере 
LUT — 
- 
aa ae 
— ES 
Sees 
gg 
re 
[m m 
A 
TTS 
Cose o c 7 
aD 
— 2 
o 
en 





Note: No entry in these columns means the coefficient was ignored. 


* 
From PRINCIPLES OF NAVAL ARCHITECTURE, SNAME, 1968, page 548. 
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TABLE 3. ASSEMSSEMNT OF THE COEFFICIENTS IN 
THE Y EQUATION 


Variable coefficient 


Series 60/5.1. 

















У m-Y v 0.209 Y9 
r mX.-Yr une ув 
V Yv -0.260 YL 
y? 1/6 Yvvv 221,450 Y2 
vr* 1/2 Yvrr 2210 Y3 
«87 1/2 Хувв 

VAU Yvu | „ес: 

уда? 1/2Yvuu | 

г (Yr-m) -0.0781 YY 4 
г? 1/6 Yrrr -0.0461 Y5 
гус 1/2 Yrvv -0.0994 Y6 
r£ 1/2 Үг46 

ráu Yru m 

rau* Угиц jos Y 
ó р 1/6 Y£ 0.020. y7 
ó v 1/2 Yérr ЭРТ 

SAU Ysu mn... 

Saus 1/2 Y Suu en. 

vr Хуга 

— Yo 0.00016 Yo 
Au Nite 
(дч? Миш € 


*No entry in these columns means the coefficient was ignored. 


*From PRINCIPLES OF NAVAL ARCHITECTURE, SNAME, 1968, page 548. 
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TABLE 4. ASSESSMENT OF THE COEFFICIENT IN 
THE N EQUATION. 























V шХ,-МУ N9 
г I,-Nr 0.0202 N8 
у Nv -0.0750 М1 
y? 1/6 Nvvv -0.3850  . N2 
vr* 1/2 Nvrr -0.3060 N3 
vs 1/2 Ку% 

vou Nvu EN аи 

vau? 1/2 Nvuu 

r (Nr-mX.) -0.0569 NN4 
r^ Шог -0.1010 N5 
rv 1/2 Nrvv -1.4200 N6 
r$ ae ber! 

rou NIU ——— 

rau? 1/2 Nruu 

Š 3 № -0.0240 №7 
| 2 1/6 1% 

бу 1/2 Кӛуч шог 

D 1/2 Nórr — 

dau Nóu —— 

£a u° 1/2 N£uu 21 

vré Nvrö 

= No -0.0003 No 
ди № NN 
(Au) Ням 2 


*No entry in these columns means the coefficient was ignored. 


*From PRINCIPLES OF NAVAL ARCHITECTURE, SNAME, 1968, page 549. 
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D. RELATIONS BETWEEN FIXED AXES AND MOVING AXES: 

So. far the reference axes X,Y, and Z are fixed to the moving ship. 
In order to determine the path of the ship, it is convenient to take the 
fixed axes to the earth X Yo: 

In figure D1., the ship is at the position A. 

The YAW ANGLE Y is the angle between the U axis and the tangent iine 
with the path at point A. | 


Then the following relationships between the fixed and the moving 


axes can be written: 


U(t) Cos Y(t)-V(t)Sin Y (t) 


X (t) = 
Е Сан Сени (е)Соз Ф(Е) 
до 
Т 
У 
И 


Figure Dl. Relations Between Fixed Axes in the 


Earth and the Ship Axes. 
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E. RUDDER DEFLECTION: 
Let Y = Rudder angle 


DRATE = Rudder deflection rate. 


TDMAX = Time at maximum rudder deflection 


Then the rudder deflection in turning maneuver can be shown in 


Figure D2; E 






5 тах|- —— -- 





"ag tmax 
Figure D2. Rudder Deflection versus Time. 


Let Rudder deflection = Ф 
DRATE = Rudder deflection rate = D_ /TDMAX 


Then the function of the rudder deflection can be written: 


i AN 
d - DRATE [елсе ад) -ВАМР Се | 57.3 
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F. THE TURNING PATH OF A SHIP: 


There are 4 phases of a turn: 





TACTICAL DIAMETER 


Ist PHASE 
EXECUTE POINT 


APR. PHASE 


+ 


Figure D3. The Turning Pain оғ а 5 


Approach phase V = r 


| 
O 
|| 
© 
< 
|| 
© 
~ 
II 
O 


First phase of turn v # 0, r £ 0, v = 0, r = 0 
Second phase of turn у% 0, r + 0, v Z 0, r £ 0 
Third phase of turn у= 0, г= 0, v # 0, r # 0 
Using dimentional equation (1) and introducing the rudder deflection 
30%. The characteristics of transient phases of a turn are shown in 


Figures D6, D7, D8, and D9. 


A-B Second phase vio, v 0, r Z 0, r #0 
B-C Third phase v #0, у = 059x507 r # 0 
Now using different rudder deflections -359 „309 the turning path 
of a ship is shown in Figures D4 and D5. When using the bigger deflection 


the steady turning radius is smaller. 
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Figure D6. Angular Acceleration about Z Axis. 
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D7. Angular Velocity about Z Axis. 
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Figure D8. Angular Acceleration about Y Axis. 
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Figure D9. Angular Velocity about Y Axis. 
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VI. SPEED CONTROL IN TURNS 


A. DYNAMIC BEHAVIOR OF THE PROPULSION PLANT: 

The propulsion plant and turning path were observed. Now combine 
both of them and add feed-back of ship speed with Gain Kl in order to 
\ 


keep the ship speed constant in turns. 


The block diagram for this part can be drawn: 


4 
# 
N 
К | 
V Ven 2) + nen | PROPULSION N HULL ve 22 
t. tX. PLANT 
7 





Figure El. Feedback Control Velocity. 


Where: V*: Reference ship speed. 
V: Ship speed. 
n*: Shaft angular speed command. 
U: Velocity about X axis. 
VÀ: Velocity about Y axis 
Kl: Ship speed gain 
n: Shaft angular speed 


6 : Rudder angle 


The block diagram of this system was shown in Appendix B. 
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B. SIMULATION METHODS 
1. Determine all initial conditions in the steady state 
RUN 1. ^ 

Let propeller shaft speed command n* be constant 
n* — 130 RPM 

And no feedback of ship speed 
Kl = 0 = 

Then determine the steady state speeds (i.e., determine 

the component of initial condition ship speed on X axis 


and on Y axis (UIC, VIC). 


RUN 2. 

Now the propeller shaft speed command n*z130 rpm making 
K1 = 8,16,24,32 respectively. Then determine the steady 
State im each case (i.e., find all IC's in the steady 
State). This guarantees that the loop is in steady 
state, since the initial conditions on other variables 
might still cause a transient. 

2. Put in a step change in the ship speed with 2 kts. increase and 

delay time 10 sec. Then V* becomes: 


V* = 20. + 2.0 step (10.) 


C. SPEED CONTROL STUDIES 

Using the inital conditions of RUN 1 and 2, a rudder command was 
applied at t = 0, causing a speed decrease due to added drag. At 
t = 10 sec, the commanded speed, V* was increased by 2 knots, stepwise. 
This caused a speed increase. The effects of changing the loop gain, 


Kl, is shown in Figure EZ. 
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D. DISCUSS THE CURVES 
1. Ship resistance 
The ship speed 
a. The ship 
then the 
b. The Wake 


then W = 


Rt, Wake fraction W, Thrust fraction t 

change was small, so that: 

resistance Rt = K.V (see Fig. B4) K = constant 
form of Rt was the same as V (Fig. E3 curve 1) 
fraction was constant (Fig. Bl) 


(Fig. E3 curve 2) 


с. In this interval the slope of the thrust fraction was 


negative. 


Then when V decreases to increases and as V 


increases t decreases (Fig. E3 curve 3) 


2. о curve, torque 


coefficient Cq, thrusting coefficient 


VP = (1-W) V 
= (1-.045) V 
= -955V 
nD 
and OMEN S — 
Таа p 


The o curve is shown in Fig. E4 curve 1 


о changed from.gs53to .920 then the slope of Ct and Cq curve was constant. 


Then the Ct and Cq curve have the same form as the J curve. 


3. Propeller action torque Qp, shaft torque Qe and propeller angular 


speed п 


Qp = 21:15 


The Qp curve was shown in Fig. E5 curve 2. 


The Wf curve and n curve are shown in Figures E6 and Ё7. 
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VII. CONCLUSION 


This thesis has presented the combined model of a propulsion plant 
and hull. DSL/360 LANGUAGE was found to be powerful for simulating ship 
motion. 

In turbine powered ships, it is necessary to predict the propeller 
response to changes in power demand. 

This problem has been investigated but the studies have been con- 
centrated on the effect of dynamics of the ship. This work is concerned 
both with the dynamic of the propulsion plant and ship dynamics of the 
ship. 

It has been shown that the complex, nonlinear dynamics of propulsion 
plant and hull can be adequately represented by an all digital computer 
model. The model permits studies of the internal dynamics of the pro- 
pulsion plant, as well as the external dynamics of the hull. 

It has been shown that the model can be incorporated in a feedback 
control system and its effects included in studies of control system 


performance. 
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VIII. RECOMMENDATIONS 


The combination of propulsion plant and hull in a single computer 
model provides a realistic tool for many studies. The following are 
topics suitable for future investigation. 

a. Course keeping, station keeping and replenishment at sea. 

b. Effect of propulsion plant dynamics on maneuvering in sea 

states (or regular waves). 

c. Effect of sea state conditions on the operation of the 

propulsion plant. 

d. Discuss variation in steady state speed with increase command 

as shown in Fig. Cl - C5. Also discuss effect of governor 


loop on steady state speed, Fig. E2. 
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APPENDIX A 


Table Look-up for Ship Resistance versus Speed: 
SHiP 


RESISTANCE 
300 Rt 


250 


200 | 
150 
100 


50 


| 
| 

| 

! 

| 

| 

| 

| 

| 

| 

{ 

! ( 


ol 5 lo 15 20 25 30 у 
FUNCTION TABR(V) 


DIMENSION VT(13),RTT(13) 

КИСА 07 205,5.0,7.5,10.,12.5,15.,17.5,20.722.5,25.5,127.5,30./ 
DATA RTT/2*0.0,7000.,13000.,25000.,39000.,57000.,80000.,1103000., 
136000. , 173000. ,225000. , 280000. / 

IF(V.GT.30) GO TO 1 

DELV 2.5 

N IFIX(V/DELV) 1 

SLOPR(RTT(N 1)-RTT (N)) / (VT (N- 1) -VT (N)) 

TABR SLOPR^*(V-VT (N)) - RTT (N) 

RETURN 

l TABR 280000. 

RETURN 

END 
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